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Abstract. Supersymmetrization of a nonlinear evolution equation in which 
the bosonic equation is independent of the fermionic variable and the system 
is linear in fermionic field goes by the name B-supersymmetrization. This 
special type of supersymmetrization plays a role in superstring theory. We 
provide B-supersymmetric extension of a number of quasilinear and fully non- 
linear evolution equations and find that the supersymmetric system follows 
from the usual action principle while the bosonic and fermionic equations 
are individually non Lagrangian in the field variable. We point out that 
B-supersymmetrization can also be realized using a generalized Noetherian 
symmetry such that the resulting set of Lagrangian symmetries coincides 
with symmetries of the bosonic field equations. This observation provides a 
basis to associate the bosonic and fermionic fields with the terms of bright 
and dark solitons. The interpretation sought by us has its origin in the classic 
work of Bateman who introduced a reverse-time system with negative friction 
to bring the linear dissipative systems within the framework of variational 
principle. 
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1. Introduction 

In the theories of elementary particles supersymmetry aims at a unified 
description of fermions and bosons i.e. of matter and interaction [1]. Under- 
standably, construction of supersymmetric theories in this case will involve 
consideration of particles of different spin and of different statistics. As op- 
posed to this, supersymmetrization of integrable nonlinear evolution equa- 
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tions proceeds by considering only the space supersymmetric invariance [2]. 
More specifically, in a superspace approach to (1 + 1) dimensional systems 
the commuting space variable x is extended to a doublet (x, 9), where 9 is an 
anticommuting variable of Grassmann type such that = 0. A superfield 
F regarded as a function of x and 9 has a very simple Taylor expansion in 
terms of 9 such that 

F(x, 9) = v{x)+9u{x) . (1) 

Here u{x) and v{x) are the component fields of F{x, 9) with f (x), the super- 
partner of u{x) and the conversly. The function F{x, 9) is fermionic (bosonic) 
if u{x) is bosonic (fermionic). 

A supersymmetric system for a given evolution equation is constructed in 
a way that will make it invariant under the transformations x ^ x — rj9 and 
9 ^ 9 + r] with t], an anticommuting parameter. Since the supersymmetric 
transformation does not depend on time we have suppressed t dependence in 
writing (1). 

For any function $ of the doublet {x, 9) it is easy to see that 

5^ = r]{de - 9d,)^ . (2) 

Thus 

Q^de- 9d, (3) 

is the generator of the supersymmetric transformation. Using $ = F{x^ 9) 
from (1) in (2) we get the transformations for the component fields as 

5v — Tju and 5u — rjVx ■ (4) 

The results in (4) give the so-called supersymmetry transformation since the 
first equation relates the bosonic field to a fermionic field while the second 
one docs the opposite. In other words, the supersymmetry transformations 
regarded as transformations in superspace exhibit the fermi-bose symmetry. 
Two successive supersymmetry transformations lead to 

— rfvx and 5^u — rfux ■ (5) 

This shows that a supersymmetry transformation is a sort of square root of 
an ordinary translation. The derivative written as 

D^9d, + de (6) 

anticommute with Q and is called the covariant superderivative presumably 
because expressions written in terms of D and superfield are manifestly in- 
variant under the supersymmetry transformation (4) . It is of interest to note 
that 

D' = d,. (7) 
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The Korteweg-de Vries (KdV) equation [3] 

Ut = -Usx + ^uu^ ■ (8) 

represents the first nonhnear evolution equation that could account for soliton 
formation [4] . The KdV equation can be extended to a supersymmetric sys- 
tem by rewriting the KdV equation in terms of the superfield and covariant 
derivative. This is achieved by multiplying (8) with 6 and then supersym- 
metrizing the result. Thus we have 

Ft = -D^F + 6DFD^F (9) 

Prom (1) and (9) we obtain equations for the bosonic and fermionic fields as 

Ut = -U3x + QuUj; (8) 

and 

Vt = -V3x + QuVx . (8') 

We note that 6{6uUx) can be regarded as a fermionic part of either 3D'^{FDF) 
or QDFD^F. Therefore, in order to construct the most general possible su- 
persymmetric extension of the KdV equation, one must consider a linear 
combination of these two terms. In writing 6{6uUx) we have not taken this 
point into consideration but chosen to work with only the term GDFD^F. 
This resulted in a set of simple equations given in (8) and (8') . Clearly, the 
bosonic equation does not depend on the fermionic variable and moreover, 
the system is linear in the fermionic field. This kind of supersymmetric ex- 
tension was originally discarded as being a 'trivial' supersymmetrization [2]. 
However, it generated a lot of interest after it was reahzed that equations like 
(8') arise in the study of superstring theory [5] . A supersymmetric extension 
in which bosonic equation do not change in the presence of fermions has now 
been given the name B-supersymmetrization. 

The aim of the present work is to provide a Lagrangian realization for 
the pair of equations which arise from B-supcrsymmetrization of nonlinear 
evolution equations. Our equations of interest are the modified KdV (mKdV) 
[6], Hunter-Saxton (HS) [7] and Camassa-Holm (CH) [8] equations. These 
equations are related to the celebrated KdV equation in someway or others. 
In this context we note that the KdV and mKdV equation are quasilinear in 
the sense that their dispersive terms are hnear while the HS and CH equations 
are fully nonlinear because their dispersive terms are nonlinear. Besides 
Hunter-Saxton and Camassa-Holm equations there exist another class of fully 
nonlinear evolution equations (FNE) introduced by Roscnau and Hyman 
(RH) [9]. These equation is apparently not related to the KdV equation. 
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In section 2 we outline how the mKdV, Hunter-Saxton and Camassa-Holm 
equations are related to the KdV equation and study their supersymmetric 
structure. We also try to present a similar supersymmetric generalization of 
the third-order RH equation. In section 3 we contemplate deriving all these 
supersymmetric equations from the action principle. 

A single evolution equation is never an Euler-Lagrange expression. We 
need to couple a given equation with an associated one to construct a phys- 
ically complete system in the sense of the variational principle. As an inter- 
esting curiosity we note that the associated equation is the fermionic partner 
of the original bosonic equation. The techanique of introducing an auxil- 
iary equation for variational formulation of physical problems has an old 
root in the classical mechanics literature. For example, in a celebrated work 
Bateman [10] noted that a dissipative system is physically incomplete and 
one needs an additional equation when one attempts to derive the defining 
equations from an action principle. In section 4 we make some concluding 
remarks. In particular, we judiciously exploit the similarity between super- 
symmetric and Bateman's dual systems to provide a physical realization for 
the bosonic and fermionic fields in terms of bright and dark solitons. 

2. Evolution equations and their supersymmetric part- 
ners Introduction 

In this section we collect a number of well-known integrable nonlinear 
evolution equations which are related to the KdV equation and study their 
supersymmetric structure. We also present a similar treatment for the non 
integrable third-order RH equation. 

(i) Modified KdV equation: The nonlinear transformation of Miura or the 
so-called Miura transformation [11] 

u^^.+e , c^ax,t) (10) 

converts the KdV equation in (8) into a modified KdV (mKdV) equation 
For breviety, introducing ^ — uwe write the mKdV equation in the form 

Ut = -Uzx + ^U^Ux . (11) 

This equation differs from the KdV equation only because of its cubic non- 
linearity. It has many applicative relevance. For example, mKdV equation 
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has been used to describe acoustic waves in anharmonic lattices and Alfven 
waves in coUisionless plasma. 

(ii)Hunter-Saxton equation: Consider the hereditary recursion operator 

R^cd'^ + X{dud-^ + u) , (12) 

for the KdV hierarchy. Here c and A are an arbitrary constants. This operator 
generates a hierarchy of integrable equations in which the first member is 

du 

— = Ru^ = cu^x^ + 3XuUx . (13) 

For c = — 1 and A = 2 we have the famous KdV equation. A second recursion 
operator can be extracted from (12) by shifting the function u to it+7, where 
7 is a constant. We thus have 

R{u + 7) = (ci9' + X{dud-^ + u)) + {C2d^ + 2A7) 

= Ri + R2 (say) . (14) 

In writing (14) we have used c = Ci + C2. It appears that the recursion 
operator R — RiR2^ generates new hierarchy of integrable equations 

ut = {RiR^^ru, . (15) 

Assuming that u — i?2^ = C2$.xx + 2A7^ the first member of the hierarchy is 

2A7^t + C2^^^t = Ci^3x + Ac2C6x + 2Ac2^a;6a; + 6A^7^^^ . (16) 

For 7 = 0, ci = 0, A = — 1 and ^ — u (16) gives the Hunter-Saxton equation 
[7] 

Uxxt + '2-UxU2x + UU^:, = . (17) 



(iii) Camassa-Holm equation: The choice ci = 0, C2 = 1, 7 = |, A = — 1 
and ^ — u (16) leads to the well known Camassa-Holm equation [8] 

l«t - «xxt + 3wWj; - 2M3;?/2a; - MMSx = . (18) 

The CH equation describes the unidirectional propagation of shallow water 
waves over a flat bottom. It can also represent the geodesic flow on the Bott- 
Virasoro group. 
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(iv) Rosenau-Hyman equation: The KdV and mKdV equations are quasi- 
linear. Here the dispersion produced is compensated by nonhnear effects 
resulting in the formation of exponentially localized solitons. As opposed to 
KdV and mKdV equations, the HS and CH equations are fully nonlinear. 
These equations do not have exponentially localized soliton solutions. In- 
stead they support peaked- and cusp-like solutions often called peakon and 
cuspon. Both HS and CH are integrable. Unlike HS and CH equations, the 
family of FNE equations proposed by Rosenau and Hyman [9] are noninte- 
grable but have solitory wave solutions with compact support. That is, they 
vanish identically outside the finite range. These solutions were given the 
name compactons. The compactons are robust within their range of exis- 
tence. In contrast with the interaction of solitons supported by KdV and 
mKdV equations, the point at which two compactons collide is marked by 
the birth of a low-amplitude compacton-anticompacton pair. The third-order 
RH equation is given by 

Ut + 3u^Ux + ()UxU2x + 'iuu^x — . (19) 

Although (11), (17), (18) and (19) appear to be structurally different from the 
KdV equation, each of them can be supersymmetrized in a rather straight- 
forward manner. In Table 1 we present results for equations of the superfield 
and supersymmetric partners resulting from B-supersymmetrization of these 
equations. For completeness we also include results for the KdV equation. 



Table 1: Nonlinear evolution equations and their supersymmetric partners 



Evolution equations 


Equations for Superfield 


Supersymmetric partners 


KdV: 
Ut + u^x - Quux = 

mKdV: 
Ut + - Gu'^Ux = 
Hunter-Saxton: 

Uxxt + UU3x + '^UxU2x = 

Camassa-Holm: 

Ut - Uxxt + ^UUx - UU^x 
-2UxU2x = 

FNE: 

Ut + Su^Ux + 6UxU2x + 
2uU3x = 


Ft + D^F - GDFD'^F = 

Ft + D^F - 6{DFfD'^F = 

D^Ft + DFD^F + D^FD'^F 
+D'^FD^F = 

Ft - D'^Ft + SDFD^F - DFD^F 
-D^FD'^F - D^FD^F = 

Ft + 3{DFfD'^F + SD^FD'^F 
+3D'^FD^F + 2DFD^F = 


Vt + V3x - QuVx = 

Vt + VSx - ^U^Vx = 

Vxxt + UVSx + UxV2x 
+U2xVx = 

Vt - Vxxt + Suvx - UV:^x 
-UxV2x - U2xVx = 

Vt + 3u'^Vx + 3UxV2x 
+3U2xVx + 2uV3x = 



Looking closely into this table we see that the superfield equations of 
KdV, mKdV and FNE equations contain only Ft as the time derivative parts 
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of these equations. In contrast to this, the corresponding time derivative 
parts of HS and CH equations are D'^Ft and Ft — D'^Ff. Understandably, the 
presence of D^Ft is associated with the mixed derivative term u^xt present 
in these equations. Note that, as expected, all fermionic fields satisfy linear 
equations. 



3. Fermionic equations from an action principle 

Our current understanding of all prototypical physical theories is largely 
based on the action principle as enunciated by Hamilton during mid 1930's. 
In the context of field theory, classical or quantum, the dynamical equations 
are obtained using the machineries of the Hamilton's variational principle. 
Thus it is expected that the bosonic and fermionic equations as given in 
Table 1 will follow from judicious use of the action principle. To achieve this 
we proceed by noting the following. 

A single evolution equation does not follow from an action principle [12]. 
When written in terms of the Casimir potential an evolution equation can 
either be Lagrangian or nonLagrangian. The example of a nonLagrangian 
evolution equation even in potential space is provided by Rosenau-Hyman 
equation [13]. Most of nonlinear evolution equation have at least one con- 
served density such that we can write these equations as 

..Mi.o. (20) 

Clearly, the KdV, mKdV and FNE equations are of the form (20). The HS 
and CH equations can be recast in a similar form by introducing n — u^x 
and m = u — u^x- 

We can make use of an elementary lemma to get a Lagrangian represen- 
tation of (20). 



Lemma 1. There exists a prolongation of (20) into another equation 

5u 



vt + ^{p[u]vx) = (21) 



with the variational derivative 

S _ ^ -iNfc ^ _ d^u 

such that the system of equations follows from the action principle 

die dxdt = . (23) 
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Here L stands for the Lagrangian density. 



Proof. For a direct proof of the lemma let us introduce L in the form 

1^ = ^(^^t - ^^t) - pH^x ■ (24) 
Prom (23) and (24) we obtain the Euler-Lagrange equations 

and 

Using (24) in (25) and (26) we obtain (20) and (21) respectively. 

Lemma 1 gives the fermionic equations for the KdV and mKdV equation 
in a rather straightforward manner. Since the HS and CH equations involve 
mixed derivative like Uxxti one needs to make use of a simple variant of 
Lemma 1 to construct the corresponding fermionic equations. To achieve 
this we write the CH equation as 

with the associated equation 

rjt + -^{P[u]vx) = (28) 

where rj — v — Vxx- The system of equations will follow from a Lagrangian 
density 

£ = -{vmt - mvt) - p[u]vx (29) 

via the Euler-Lagrange equation 

d (P . dC d dC 5C 

Jtd^^d^t JtWt '^^ 

obtained from (23) and (29). Here (j) is either u or v. A similar treatment 
also applies for the HS equation. 

Application of Lemma 1 to the FNE equation gives the superpartner 

Vt + 3u^Vx + 2uv3x = . (31) 
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This equation does not agree with the result given in Table 1 presumably 
because (19) is nonLagrangian even in the potential representation. In the 
recent past one of us (BT) [13] found that the Helmholtz solution of the 
inverse problem [14] can be judiciously exploited to construct Lagrangian 
systems from (19), which support compacton solutions. For example, the 
equation 

Ut + Su^Ux + 8uxU2x + 4:UU3x — (32) 

when written in terms of Casimir potential w{x, t) = u{y, t)dy was found 
to result from an action principle and support compacton solution. The 
fermionic equation corresponding to (32) can be obtained as 

Vt + M^Vx + ^UxV2x + ^U2xVx + ^UV^x = . (33) 

It is easy to verify that (33) can also be obtained by the use of Lemma 1. Thus 
we infer that all fermionic equations resulting from B-supersymmetrization 
of Lagrangian system of nonlinear evolution equations can be obtained by 
using our Lemma 1. 

About a decade ago Kaup and Malomcd [15] sought an application of the 
variational principle to nonlinear field equations involving dissipative terms. 
These authors demanded that the Lagrangian density 

C — v{x,t) X (equationfor it(x, t)) . (34) 

should be minimize the action functional. An important consequence of 
writing (34) is that in the presence of the auxiliary field v{x^ t) the resulting 

set up Lagrangian symmetries coincides with the symmetries of the evolution 
equation for u{x,t). In point mechanics this type of generalization for the 
traditional Noetherian symmetry was introduced by Hojman [16]. 
For the KdV equation (34) reads 

C^v{ut + 6uux + U3x) ■ (35) 

Clearly, for C in (35) the Euler-Lagrange equation in v{x,t) gives the KdV 
equation and more significantly we get the corresponding fermionic equation 
from the Eulcr- Lagrange equation for u{x,t). Starting from (34) we can 
verify the results given in Table 1 for the mKdV, HS and CH equations. 
Consistently with our previous analysis the use of (34) does not reproduced 
the tabulated result for the RH equation. Rather we find that (33) is the 
supersymmetric partner of (32) which is Lagrangian in the potential space. 
Thus there must exist a relation between the Lagrangian densities in (24) 
and (34), For the KdV field the gauge terms ^ (l"^"^) + ^ {3u'^v) + ^ {u2xv) 
can be added to (24) to get (35). A similar conclusion also holds good for 
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other Lagrangian equations in Table 1 implying that the Lagrangians in (24) 
and (34) are gauge equivalent. 



4. Conclusion 

Nonlinear evolution equations do not admit Lagrangian representation in 
terms of the field variables because the Frechet derivative Dp of any equation 
written as ut = P[u] is non self-adjoint such that the Helmholtz theorem [14] 
becomes inapplicable. Similar nonLagrangian systems also occur in point 
mechanics. For example, it is well known that there is no direct method of 
applying variational principles to nonconservative systems, which are char- 
acterized by friction or other dissipative processes. In fact, one can not write 
a time-independent Lagrangian even for a simple system like the damped 
Harmonic oscillator represented by 

iji + Xqi + = (36) 

since it is non self-adjoint. Here A is the frictional coefficient and cu, the 
natural frequency of the oscillator. Batcman [10] showed that while (36) is 
non Lagrangian, the variational principle can be applied to a dual system 
consisting of (36) and the equation 

q2 - Xq2 + (^'^q2 = . (37) 

Equations (36) and (37) can be obtained by using the Lagrangian 

L = q2{qi + \qi + uj'^qi) (38) 

in the Euler-Lagrange equations 



dL_ 
dq2 



and 

dL d ( dL 



(39) 



= . (40) 




dqi dt \dqi 

Note that L in (38) can be converted to the well-known Morse-Feshbach [17] 
Lagrangian by adding appropriate gauge terms. 

Equation (36) represents the motion of a simple pendulum embedded in a 
fluid which opposes the motion through frictional forces proportional to the 
velocity. As a result the energy is drained out of the system. The complemen- 
tary equation (37) represents a second physical system which absorbs energy 
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dissipated in the first such that (36) and (37) taken together can be viewed 
as a conservative system to have an exphcitly time independent Lagrangian 
representation. Thus the physical interpretation of Bateman dual system is 
fairly simple and straightforword. This is, however, not the case with the 
pair of equations obtained by B-supersymmetrization and finally shown to 
form a Lagrangian system. Nonintegrability of these equations [2,5] might 
be one the reasons for that. 

The single soliton solution of the KdV equation in (8) is given hj u = 
—2k'^ sech? k{x — ikH). The solution has amplitude —2k'^ and moves to 
the right with speed Ak"^. For this value of u we have found v in (8') as 
V = tanh? k{x — Ak'^t). The vave represented by v is of constant amplitude 
and moves with same velocity as that of u. Moreover, \u\ with a peak while 
1^;! is a notch. Looking from this point of view u and v may be regarded 
as bright and dark soliton respectively. The bosonic and fermionic fields of 
(32) and (33) are associated with compaction and anticompacton expressed 
in terms of appropriate Jacobi elliptic function. This kind of realization 
for bosonic and fermionic fields marks a point of departure from our under- 
standing of quantum field theory in that the particles and antiparticles are 
not superpartner. 
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